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Abstract: 

I ^ . We study the diffusion-limited process A + A ^ A in one dimension, with finite 

QQ . reaction rates. We develop an approximation scheme based on the method of Inter-Particle 

Distribution Functions (IPDF) , which was formerly used for the exact solution of the same 
>■ ■ process with infinite reaction rate. The approximation becomes exact in the very early 

time regime (or the reaction-controlled limit) and in the long time (diffusion-controlled) 
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O ■ asymptotic limit. For the intermediate time regime, we obtain a simple interpolative 

Q . behavior between these two limits. We also study the coalescence process (with finite 

^ . reaction rates) with the back reaction A ^ A + A, and in the presence of particle input. In 

each of these cases the system reaches a non-trivial steady state with a finite concentration 
of particles. Theoretical predictions for the concentration time dependence and for the 
IPDF are compared to computer simulations. 
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1. Introduction 

In recent years much effort has been dedicated to Diffusion-limited reactions in low 
dimensions. ^~^ Most research has focused on the bimolecular reaction A + B ^ inert, ^ 
and on one-component coalescence, A + A ^ A,^~^ and annihilation, A + A ^ 0.^'^ 
The last two systems were solved exactly in one dimension and are especially useful in 
elucidating the anomalous kinetics of diffusion-limited reaction processes. Naturally, simple 
generalizations and extensions of these basic processes are of much interest, because of 
the possibility that they may also be solved exactly. Indeed, the coalescence process, 
A + A -^ A, has been also solved together with the back reaction, A ^ A + A, and under 
input of A particles,^ as well as for systems of finite size^ and with inhomogeneous initial 
conditions.^ 

In the diffusion-limited coalescence process reactions occur at an infinite rate, taking 
place immediately upon the encounter of any two particles. An obvious generalization 
would be to make the reaction rate finite. This more physical process would also ex- 
hibit an interesting crossover from a classical, reaction-limited behavior, at early times, 
to a diffusion-limited regime, in the long time asymptotic limit. Surprisingly, this sim- 
ple minded generalization makes it difficult (perhaps even impossible) to solve the model 
exactly. 

Simulations have been performed, both on the one-component annihilation process 
and on the coalescence process in one dimension, with a finite reaction rate.^'^"^^ The 
simulations show three different regimes: (a) an early time regime where the particles 
merely diffuse with a negligible change in the initial concentration, (6) an intermediate 
regime, in which the concentration of particles decays faster than the diffusion limited 
case, c ~ 1/v^, but not quite as fast as the classical limit, c ~ 1/t, and (c) a long 
time, diffusion-limited regime, where the system behaves exactly as if the reaction rate 
were infinite. Based on these simulation results, the intriguing possibility that in the 
intermediate regime the concentration decays anomalously in time, with a power dependent 
on the reaction rate, was raised. ^° Subsequent theoretical work by Privman et al}"^ and 
Hoyuelos and Martin^^ has suggested that the kinetics in the intermediate regime can be 
explained as an interpolation between the classical and the diffusion-controlled limits. 

The analysis of Privman et al}"^ is based on the concept of inter-particle distribution 
functions (IPDF), which had been introduced earlier for the exact solution of the coa- 
lescence process with infinite reaction rate.^ In the IPDF method a diffusion equation is 
derived for the probability that two nearest particles are a distance x apart (the IPDF), 
and the reaction is represented by absorbing boundary conditions at the origin. Privman 
et al. have replaced this boundary condition with a radiative boundary condition, in or- 
der to approximate the finite reaction rate. Their approach yields an elegant, qualitative 
understanding of the three regimes discussed above. Hoyuelos and Martin^^ have found 
an interpolation formula for the challenging intermediate regime, which fits the simulation 
data reasonably well. However, the derivation is largely phenomenological and it requires 



a scaling ansatz which we find physicaUy obscure. 

In this paper, we study the one-dimensional coalescence process, A + A ^ A, with 
finite reaction rate. Our approach is based on the IPDF method but differs from that of 
Privman et al. in important details. In particular, it allows us to derive a closed analytic 
expression for the concentration decay at all times. For the intermediate regime, an inter- 
polative formula between the classical and diffusion-controlled limits emerges as a natural 
consequence of the approximation. We also study the coalescence process with back reac- 
tions, A ^ A + A, and with particle input, ^ A. In each of these two cases the system 
arrives at a non-trivial steady state where the concentration is finite. Finally, we study the 
IPDF's themselves, which show different characteristics at different time regimes and in the 
stationary cases. Our theoretical results compare quite well with computer simulations. 

The rest of the paper is organized as follows. Our model and a physical approximation 
based on the method of IPDF's are introduced in Section 2. In Section 3, we perform a 
mathematical approximation that enables us to solve the equations derived in Section 2 
in closed form. In the same section, we compare between simulations and analytic results. 
The two cases of non-trivial steady state are explored in Section 4, both through our 
analytic approach and through computer simulations. We conclude with a discussion, in 
Section 5. 

2. The reaction model and the IPDF method 

Our model is defined on a one-dimensional lattice with lattice spacing Ax. Each site 
can be either empty (o) or occupied (•) with one particle. The particles move randomly, 
independently of each other, to a nearest neighbor site with a hopping rate D /[/S.xY (to 
each side). On long length and time scales this yields normal diffusion with diffusion 
coefficient D. When a particle hops onto a site which is occupied, coalescence takes place 
with probability /c, while with probability (1 — /c) the particle is reflected back to its 
original position and no reaction takes place. Thus, the probability k controls the rate of 
the reaction 

A + A^A. (1) 

When /c = the reaction rate is zero. The particles merely diffuse, bouncing off each other 
and no reactions take place. In the other extreme, when k = 1 reactions are immediate, 
i.e., the reaction rate is infinite. This is the purely diffusion-limited case. When k is very 
small but finite, one expects to find a regime which is dominated by the slow reaction rate. 
This is the classical, reaction-limited case. 

To solve the system, we follow the IPDF method, previously used for the exact solution 
of this model in the case of /c = 1.^ We define, as usual, En{t), the probability that a 
randomly chosen segment of n consecutive sites is empty, i.e., contains no particles. The 
probability that a site is occupied is 1 — Ei = C, and the density of occupied sites is 

c(t) = (1 - Ei{t))/Ax = C{t)/Ax. (2) 



En gives the probability that, say, sites 1 through n are empty, while -En+i gives 
the probability that sites 1 through n are empty and also site n + 1 is empty. Thus, the 
probability that a segment of n sites is empty, but that the adjacent (n + l)th site is 
occupied, is Prob(o o ■ ■ ■ o •) = En — En+i- It can also be shown^ that 

En-l — '^En + En+l = Cpm (3) 

where Pn is the probability that nearest neighbor particles are exactly n sites apart (this 

is the IPDF, from which the method derives its name). In particular, the probability of 

finding two adjacent occupied sites is exactly given by Prob(««) = 1 — 2Ei + E2. 

When the reaction probability /c = 1, one can write down an exact equation for the 

evolution of E^:^ 

2D 

The first term on the r.h.s. describes the creation of an empty n-sites interval, when a 
particle at the inner edge of the interval hops out. The second term describes annihilation 
of an empty interval, when a particle at the outer edge of the interval hops in. 

When the reaction probability /c < 1, hopping out of the interval is not always pos- 
sible. The target site may be occupied, in which case hopping (and coalescence) will be 
disallowed, with probability 1 — k. To account for this effect, we require the probability of 
finding intervals of n — 1 consecutive empty sites, followed by two occupied sites. Since this 
probability cannot be expressed in terms of the En exactly, we propose the approximation 



dtEn = j^^—^{{En-l - En) " (K " K+l)}- (4) 





OProb(««) {l-2Ei + E2){En-i-Er, 



' Prob(«) 1 - El 

With this approximation, the evolution equation for arbitrary k is 



(5) 



dtEn = J^iEn-l - 2En + En+l) - -|^(1 - k) ^ ^^'^^^ {En-1 - En). (6) 

where the correction due to failed coalescence attempts is represented by the last term. 
Notice that when k = 1 this properly reduces to Eq. (4). 

Eq. (6) is valid for n > 1. For n = 1 we have the exact equation: 

^tE^ = J^k{l-2E^+E2), (7) 

which simply states that sites become empty at the same rate as particles coalesce. Com- 
paring Eqs. (6) and (7) we see that they may be combined, by requiring the boundary 
condition 

Eo{t) = 1. (8) 



Eq. (6) may thus be extended to the case of n = 1. A second boundary condition is 

En{t) =0, as n ^ oo, (9) 

since an infinitely long segment will always contain particles, as long as the concentration 
is finite. 

The only approximation made here is in Eq. (5), where some correlations between the 
probable state of consecutive intervals are neglected. We argue, however, that Eq. (6) is 
asymptotically correct both in the early and long time regimes, and hence it may provide a 
reasonable interpolation for the intermediate time regime. If the starting configuration of 
the system is random, then E^ = E^ and the state of consecutive intervals is uncorrelated. 
This situation will persist, and Eq. (5) will hold, until the concentration drop is noticeable 
{i.e., the end of the early time regime) for it is only reactions that induce correlations 
(in fact, diffusion randomizes the system). After very long times, on the other hand, the 
concentration of particles becomes very small. As a result, adjacent occupied sites become 
extremely rare and the correction term in Eq. (6) eventually becomes negligible. Eq. (6) 
then degenerates to the case of infinite reaction rate and it does no longer matter how 
imprecise the correction term is. 

3. Integration of the evolution equation 

There exist various techniques to solve Eq. (6), with the boundary conditions of 
Eqs. (8) and (9). Perhaps the most straightforward method is numerical integration, 
which requires discretization of the time variable in Eq. (6). A second approach, which 
worked well for the case of /c = 1, is passing to the continuum limit. ^ This is achieved 
by defining the spatial coordinate x = nAx. The probabilities En{t) are replaced by the 
function E{x,t). Letting Ax — > 0, Eq. (6) is replaced by 

dtE = 2DdlE - 2D{1 - k)^^^^^d^E. (10) 

(JxE\x=0 

and the boundary conditions of Eqs. (8) and (9) become 

£;(0, t) = 1, and E{x ^ oo, t) = 0. (11) 

We see that the finite reaction rate gives rise to a non-local, non- linear term in Eq. (10). 
However, notice that 9^£'|a;=o/(^x-E'|x=o = '^(t) is a function of time only. Thus, in principle 
one can proceed by Laplace-transforming Eq. (10) with respect to the spatial variable x, 
and then determine uj{t) in some self-consistent way. Unfortunately, this procedure leads 
to complicated expressions and one is forced to resort to series expansions, limiting the 
solution to a few asymptotic results. 

Here we propose an alternative approach, based on an approximation of Eq. (6). We 
emphasize that this approximation is a mere mathematical convenience, designed to enable 



us to obtain a solution to Eq. (6) in closed form. The physical approximation made in 
Eq. (5) is the real focus of this paper. Indeed, Eq. (6) can be solved to any degree of 
accuracy employing numerical methods. 

We first sum Eq. (6) over the index n, from 1 to oo, to yield 

n=l ^ ' ^ ' 

where we have used the boundary conditions of Eqs. (8) and (9). The r.h.s. can be made 
a function of only Ei, with the help of Eq. (7). For the l.h.s. we make the approximation 
^ En ~ A/(1 — Ei), where A is a constant. The motivation for this is that in the long time 
asymptotic limit E^ ~ 1 for all n up to a characteristic (n) = 1/(1 — Ei), and falls sharply 
to zero for n > (n) . More precisely, in the long time asymptotic limit the reaction proceeds 
as if k is effectively 1, in which case we know that A = 2/7r, exactly.^ The approximation 
lies in the fact that we assume A to be constant at all times. Indeed, the variation in A is 
quite small; at the beginning of the process, when the distribution is random, A = 1. Let 
us then assume A = Ijix (to match the exact long time asymptotic solution) to hold true 
throughout the process. Eq. (12) then becomes 

" 2 c+i^^. (13) 



dr TiC kC dr 

where r = (2Z)/(Ax)^)t is a dimensionless time variable, and we have used 1 — i?i = C 
(Eq. 2). The solution to Eq. (13) is 



C= ^ , , . —. (14) 

where Co = C{t = Q). 

In Fig. 1, we plot the concentration decay as computed from a numerical integration 
of Eq. (6), from the analytical expression of Eq. (14), and from computer simulations, for 
the same choice of parameters {D and Co) and for various choices of k. The agreement 
between Eq. (14) and the numerical integration is excellent — little is lost in the 'mathe- 
matical' approximation. More importantly, there is good agreement between theory and 
simulations: the early and late time regimes match almost perfectly, and in spite of differ- 
ences of up to 9 % in the intermediate regime, the slope of the curves and the crossover 
times are almost identical. 

From Eq. (14), we see that the intermediate time regime is merely an interpolation 
between a classical decay, C ~ 1/t, and a diffusion-limited decay, C ~ l/v^. We can use 
our result to estimate the crossover times. We first expand C{t) in powers of t: 
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We obtain the crossover time between the early time regime and the intermediate regime, 
Ti, by requiring that the term hnear in r be a finite fraction, e, of Cq: 



2 1-k, 

'ttCJ^ kCo 



Ti=e{-^ + ^:^). (16) 



Even in the diffusion-hmited case, when k ^ 1, there is an early time regime where reactions 
go unnoticed. Since coalescence is immediate, the crossover time equals the typical time 
that two nearest particles will take to reach each other. The average distance between 
particles is I/cq, and since the particles diffuse, ti ~ 1/Dcq,^^ in agreement with the 
estimate above. The second term on the r.h.s. predicts that the crossover time ri will 
increase proportionally to I/Zccq, which is characteristic of the classical limit (classically, 
dc/dt = —kc^ and c = co/(l + kcot)). 

Next, we expand Eq. (14) in powers of l/^/r? 

\/7lT 2KT T^! 



Notice that the leading term corresponds to the long time asymptotic limit, where the 
concentration decays as c = l/yhxDt (as in the exact solution for k = 1). Interestingly, 
the second term, too, does not retain any memory of the initial density, Co, but has some 
/c-dependence. Comparing the leading term to the first correction we get an estimate for 
r2, the crossover time between the intermediate regime and the long time regime: 

^2 - -^^^- (18) 

This can be explained heuristically, as follows. The long time asymptotic regime occurs 
because the reaction probability k effectively renormalizes to 1: If the density of particles 
is C, the number of sites between neighboring particles is on the average 1/C. Because 
the particles diffuse, it takes them of the order of 1/C^ steps to meet each other. During 
this time, each of the 1/C sites is visited of the order of (1/C^)/(1/C) = 1/C times. In 
particular, two neighboring particles will collide about 1/C times before wandering off each 
other to interact with other particles. This means that two neighboring particles will react 
almost surely, before meeting other partners, if (l/C)/c ~ 1. That is, the crossover will 
occur when C ~ /c, or, since in the long time asymptotic regime C = l/y^7rr, T2 ~ I/ttA;^, 
in agreement with Eq. (18). 

To obtain the IPDF, we integrate Eq. (6) numerically and then use Eq. (3). At time 
t = 0, we start with a random distribution of particles, so that Pnifi) = (1 ~ Cq)"^ /Cq. 
In the continuum approximation, this can be written in the scaling form p(^) = exp(— ^), 
where ^ = c{t)x. As the reaction proceeds, the likelihood to find nearby particles decreases, 
due to the coalescence process. Thus, coalescence gives rise to an effective repulsion. In the 
long time asymptotic limit, when k effectively renormalizes to 1, the IPDF arrives at the 



stationary scaling fornip(^) = (7r^/2) exp(— 7r^^/4).^ During the intermediate time regime, 
the IPDF makes a smooth transition between these two hniits. In Fig. (2), we plot p{^) 
as obtained from numerical integration and compared to computer simulations at various 
stages of the process. As might be expected, the agreement is worst in the intermediate 
time regime. 

4. Coalescence with non-trivial steady states 

Until now, we have treated the case where in the long time asymptotic limit the 
concentration drops to zero. We now want to discuss situations with a non-trivial steady 
state. It is important to examine such situations because they impose a stricter test on 
the validity of our approximation. 

Back reactions 

We consider first the case where the reaction (1) is reversible:^ 

A + A^A. (19) 

In the reverse process, A ^ A + A, a, particle gives birth to another at an adjacent 
site, at rate v/Ax (this means rate vjlAx on either side of the original particle). The 
corresponding evolution equation is 

^'^- = (A^^^--^ ' ^^^ " ^"+^^ ' (A^^^ " ^^ \-E, '^^-' - ^-^ (20) 

V 
— -^{En — En+l), 

where the last term describes the annihilation of an n-sites empty interval due to a birth 
event from a particle at its outer edge. 

In the steady state, the l.h.s. is equal to zero and Eq. (20) becomes a recursion relation 
for the En- The solution can be found by assuming a steady state of maximum entropy: 
since the process is reversible, the steady state is an equilibrium state and hence the IPDF 
is a Poisson distribution, E^ = i?f . Taking into account the boundary conditions of 
Eqs. (8) and (9), we find 

^- = ^2Dk + vAx^^- ^^^^ 

and 

_1- El _ V 

^" ~ Ax ~ 2Dk + vAx' ^^^^ 

Notice that although Eq. (20) contains an approximation, the corresponding steady-state 
equation is exact: when the IPDF is completely random the approximation of Eq. (5) 
becomes exact. Thus, the result of Eqs. (21) and (22) is exact. This is well confirmed by 
simulations. 



Particles input 

Consider now the case of a random, steady input of particles. At each time step empty 
sites become occupied with probabiUty RAx. That is, R is the increase in concentration 
per unit time, due to input. ^ The evolution equation is 

dtE^ =^^{E^_^ - 2E^ - £;„+i) - -|^(1 - k) ^ ~f \+ ^^ (£;„-i - E^) 

[AxY (Ax)^ 1 - El (23) 



- RnAxE; 



ni 



where the last term represents the annihilation of an n-sites empty interval due to input. 
Here the steady state limit is less simple than for back reactions. Although the in- 
put sustains a steady concentration of particles at the long time asymptotic limit, this 
stationary state is not a true equilibrium state (the input process is not the reverse of coa- 
lescence) and the particles are more ordered than in a Poisson's distribution. Nevertheless, 
Eq. (23) with dtE^ = still yields an approximate recursion relation for the stationary 
En. This can be solved exactly, but the solution is rather cumbersome. Instead, it is more 
enlightening to consider the continuum limit of the steady state equation: 

= 2DdlE - 2D(1 - k)ud^E - xRE, (24) 

where to = d'^E\x=Q / dxE\x=o = —dlE\x=o/cs is now a constant. To determine u, let us 
look at the discrete steady state equation for n = 1, 

2D 

0=-—-:rk{l-2Ei + E2)-AxREi. (25) 

It simply equates the rates of input events and coalescence events in the steady state. The 
continuum limit of Eq. (25) is too drastic in that it yields zero for each of these rates (and 
also a; = 0). 

A somewhat inelegant, but effective way around this is to retain AxR finite, so that 
uj = —RAx/2DkCs ^ 0. Then, the solution to Eq. (25) with the boundary conditions of 
Eq. (11) is 

E{x) = exp ..r _2/3/. Wi : 26 

KCs Ai[r ^/"^/{kCs) \ 

where k = 4Dk/{l — k)RAx and r = R/2D. From the relation Cg = —dxE\j;^o, we then 
obtain a transcendental equation for the steady state concentration, Cg: 

, - 1 ^,/3 Ai-[r-^/3/(.c.)^] 
' KCs Ai[r-2/3/(Kc,)2] ■ ^ ' 

In Fig. 3, we plot Cg as a function of the reaction probability /c, for fixed r (r^' ^ = 0.04), 
as obtained from computer simulations and from numerical integration of the discrete 
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steady state equation. The agreement between simulations and theory is quite good. Also, 
for the range shown, the agreement between Eq. (27) and the numerical integration is 
better than 4 %. 

For small concentrations, Eq. (27) itself may be simplified. If c <^ t^'^/k, then 
Ai'[r-2/3/(KC«)2]/Ai[r-2/3/(KC5)2] ~ Ai'(0)/Ai(0) = 0.72901 . . ., and the equation reduces 
to a simple quadratic, with the solution 



1/. ,,/,,, 2(l-fc)A^ N - _ Ai'(O) ,3 



Here Cg is the steady state concentration when /c = 1, and is exact. ^ In the range shown 
in Fig. 3, Eq. (28) agrees with the Cs obtained from numerical integration of Eq. (23) 
to within 10 %. We emphasize, again, that the mathematical approximations made here 
are not essential to our method, but are merely done to obtain simple final expressions, 
to better understand the consequences of the physical approximation made in writing 
Eq. (23). 

Finally, the IPDF itself is obtained from p[x) = {l/cs)d'^E/dx^, and using Eq. (26) 
with the Cs found from numerical integration (or any of the approximation formulae). In 
Fig. 4, we compare between the IPDF's obtained in this way and from computer simulations 
for r^/^ = 0.04 fixed, and various values of k. The tail of p{x) falls off as exp(— x^/^) and 
for moderately large k p{x) shows a maximum. Thus, the system is more ordered than in 
the case of back reactions, but less ordered than in the case of pure coalescence, where the 
tail of the IPDF decays as exp(— x^). The agreement between theory and simulations is 
best for large k. From Figs. 3 and 4 one can see that the IPDF is a much more sensitive 
test for approximations than the concentration alone. 

An interesting result is that p{x — 0) is zero only when k = 1. This is in contrast to 
the case of pure coalescence, where in the long time asymptotic limit p{0) = 0, regardless 
of the value of k. Recall that in the latter case k effectively renormalizes to 1 when 
the concentration decreases. This cannot happen in the case of back reactions, when c 
approaches a stationary finite value. Then, if the reaction rate is small enough it can 
overcome the strong effective repulsion between nearest particles. 

5. Summary and discussion 

We have introduced an approximation, based on the IPDF method, that allows us 
to draw analytically simple descriptions of the coalescence process A + A ^ A in one 
dimension, with finite reaction rates. The intermediate time regime is convincingly shown 
to be a crossover behavior between the classical, reaction-controlled limit, in which the 
concentration decays as c ~ 1/t, and the diffusion-controlled limit, where c ~ 1/Vt. This 
result, summarized in Eq. (14), is similar to that of Hoyuelos and Martin. ^^ It requires 
no fitting parameters and no approximations other than the standard approximation of 
Eq. (5). 
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The approximation yields exact results for the equilibrium state, reached in the pres- 
ence of the reverse reaction A ^ A + A. This is because in the equilibrium state the 
fundamental approximation of Eq. (5) becomes exact. Surprisingly, we obtain rather accu- 
rate results for the stationary state of the process with input, in spite of the strong spatial 
correlations that do evolve. 

The present approach is similar to that of Privman et a/.^^ Their equation for the 
IPDF, 

dtpix, t) = 2Ddlp{x, t) + 2D(1 - /c)p(0, t)dp{x, t), (29) 

is closely related to our Eq. (10), through p{x,t) = d^E{x,t)/c{t). The difference is that 
Prob(oo- ■ -o**) ^ Prob(oo- ■ ■o«)Prob(«) is used, instead of our Eq. (5). Although Eq. (5) 
is a bit more sensitive to short range correlations, we believe that this difference is trivial. 
The advantage of our approach is more likely in that in writing an evolution equation for 
E{x,t), we are able to use the exact boundary condition E{0,t) = 1. In contrast, there is 
no such constraint on p{x = 0, t) of Eq. (29). 

Regarding reverse reactions and input of particles, we have discussed here only the 
stationary limits and have not addressed the time-dependent relaxation to the steady-state. 
A straightforward separation of variables and decomposition into eigenvalue equations is 
not possible, because of the non-linear nature of the correction term arising from the 
finite reaction-rate. In the case of back reactions and infinite reaction-rate, there is a 
dynamical phase transition in the characteristic relaxation time.^ It will be interesting to 
see how is this transition affected by finite reaction rates. Is there a range of reaction- 
rates for which the transition disappears? The approximation method used here has been 
successfully employed, with slight modifications, for the study of the diffusion-limited 
processes 3 A -^ 2 A and 3A — » A,^^ and for the contact process in one dimension. ^^ It may 
be worthwhile to further explore its potential. For example, Eq. (5) can be systematically 
improved, by writing an hierarchy of evolution equations for the state probabilities of 
finite size intervals and introducing a truncation scheme at a later stage. These intriguing 
questions are left for future work. 
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Figure Captions 

Figure 1: Concentration decay as a function of time for tiie coalescence process as obtained 
from computer simulations (solid line), and from Eq. (14) (broken line). The different 
curves represent different values of the reaction probability; k = 0.005, 0.01, 0.02, 0.04, 
0.08, and 0.16 (from top to bottom). Shown also are results from numerical integration of 
Eq. (6) for k = 0.02 (circles). 

Figure 2: The scaled IPDF, p{cx), for the coalescence process for k = 0.02, as obtained 
from computer simulations (circles) and from numerical integration of Eq. (6) (solid line). 
Results are shown for t = 0, 100, 1000, and 10000 (p(0) decreases with time) . The 
crossover times for this particular process are roughly at ti = 15 and ^2 — 4000. Shown 
also is the IPDF in the long time asymptotic limit (broken line). 

Figure 3: The stationary concentration, c^, as a function of the reaction probability, k 
for coalescence with input of particles, with r = (0.04)"^. Numerical integration of Eq. (23) 
(circles) is compared to computer simulation results (solid line). 

Figure 4: The scaled IPDF, p{csx) for the same process as in Fig. 3 for k — 0.02, 0.04, 
0.08, 0.20, and 0.40 (p(0) is smaller for larger fc), as computed from the second derivative 
of Eq. (26) (broken line) and compared to computer simulations (solid line). 
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